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In this paper, we study the geometry of reduced density matrices for states with symmetry-protected topo¬ 
logical (SPT) order. We observe ruled surface structures on the boundary of the convex set of low dimension 
projections of the reduced density matrices. In order to signal the SPT order using ruled surfaces, it is important 
that we add a symmetry-breaking term to the boundary of the system—no ruled surface emerges in systems 
without boundary or when we add a symmetry-breaking term representing a thermodynamic quantity. Although 
the ruled surfaces only appear in the thermodynamic limit where the ground-state degeneracy is exact, we an¬ 
alyze the precision of our numerical algorithm and show that a finite system calculation suffices to reveal the 
ruled surface structures. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Dd, 03.67.Mn 

I. INTRODUCTION 

The idea of reading physical properties from the geometry 
of the underlying convex bodies arising naturally from quan¬ 
tum many-body physics has been examined and re-examined 
many times from different perspectives in the literature [1- 
6] . It appeared in the context of the W-representability prob¬ 
lem in quantum chemistry and, more generally, the quantum 
marginal problem in quantum information theory [1, 7]. Re¬ 
cently, the approach received revived interests and was used 
in the investigations of quantum phases from the convex ge¬ 
ometry of reduced density matrices or their low dimension 
projections [3-6, 8]. 

Consider a many-body local Hamiltonian H = ITj 
where each term Hj acts non-trivially on a set of particles Sj . 

Usually, each Sj contains only a constant number of neigh¬ 
boring particles and defines the geometry of the local interac¬ 
tions of the system. The convex set (a convex set in Euclidean 
space is the region that, for every pair of points within the re¬ 
gion, every point on the straight line segment that joins the 
pair of points is also within the region. For instance, a solid 
sphere forms a convex set) of quantum marginals, the collec¬ 
tions of reduced density matrices (pj), is of fundamental im¬ 
portance to quantum many-body physics. Here, each pj in 
the same tuple is the reduced state on Sj of a common many- 
body state. The convex set of quantum marginals is indepen¬ 
dent of the particular form of the Hamiltonian H, but depends 
only on the geometric locality of the system defined by 5'j’s. 

Should we have a complete characterization of the quantum 
marginals, most problems of many-body physics would be¬ 
come extremely easy. As expected, however, the structure of 
the set of quantum marginals is rather complicated. 

To study quantum phase transitions, it seems appropri¬ 
ate to investigate a further coarse-grained set of the quan¬ 
tum marginals. The usual situation one often faces when 


considering quantum phase transitions at zero temperature is 
the following: the Hamiltonian of the system has the form 
H = JiHi J2H2 and one is concerned with the change 
of the properties of lower energy states of H as Ji and J2 
change. The terms Hi and H 2 now act on a large number of 
particles, but they are still sums of local terms. Let the convex 
set0({i7j}) C be the set of all points (tr(i7ip), tr(i72p)) 
for p ranging in the set of all possible many-body states. This 
set is a low dimension projection of the convex set of quantum 
marginals and is hopefully much easier to analyze. 

It is obvious that for any ( 0 : 1 , 0 : 2 ) in Q{{Hj}), Jiai > 
Eo{H), the ground state energy of H. This means that the 
Hamiltonian H can be thought of as the supporting hyper¬ 
plane of 0({i7j }), and the change of parameters Ji, J 2 can be 
visualized as the change of the supporting hyperplane, mov¬ 
ing around the convex set. The intersection of this hyperplane 
with &{{Hj}) corresponds to the image of the ground state 
of H on the boundary of 0({i7j}). A flat portion on the 
boundary of Q{{Hj}) signatures the flrst-order phase tran¬ 
sition. However, for continuous phase transitions, the geom¬ 
etry of &{{Hj}) alone does not convey any informative sig¬ 
nals [8]. 

Recently, the convex geometry approach was employed in 
the study of quantum symmetry-breaking phases [9]. For 
a Hamiltonian H = JiHi J 2 H 2 with certain symmetry, 
one adds a third, symmetry-breaking, term to the Hamil¬ 
tonian and consider H = JiHi f- J 2 H 2 + JsH^. The 
authors plotted the convex set 0({i7j}) C of points 
(tr(i7ip),tr(i72p),ti*(^3P)) analyzed the geometry of 
its boundary. On this set, the emergence of ruled surfaces 
on the boundary is observed (a ruled surface is a surface that 
can be swept out by moving a line in space, or equivalently, 
for any point on the ruled surface there exists a line passing 
through this point that is also on the surface. For example, 
the curved boundary of a cylinder is a ruled surface). The au- 
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thors of [9] argued that the existence of those ruled surfaces is 
a defining property of symmetry-breaking and can be used to 
signal symmetry-breaking phase transition. 

Interestingly, the observation that ruled surfaces on the 
boundary of certain convex body can explain phase transitions 
dates back to Gibbs in the 1870’s [10-13], even though the 
convex bodies under consideration in classical thermodynam¬ 
ics and quantum many-body physics are rather different. It 
indicates that the convex geometry approach is a rather funda¬ 
mental and universal idea. 

In the present paper, we study the phenomenon of the 
emergence of ruled surface on the boundary of the convex 
set 0({i7j}) C for one dimensional (ID) symmetry- 
protected topological (SPT) ordered systems. An SPT or¬ 
dered state is a bulk-gapped short-range entangled state with 
symmetry protected nontrivial boundary excitations [14]. The 
well known two dimensional and three dimensional topologi¬ 
cal insulators [15-20] are free fermion SPT phase protected by 
time reversal and U{1) charge conservation symmetry, whose 
boundary remain gapless as long as the symmetries are pre¬ 
served. SPT phases also exist in interacting systems. A typi¬ 
cal example of interacting bosonic SPT phase is the ID spin-1 
Haldane chain [14, 21-24], whose degenerate edge states are 
protected by time reversal, or spatial inversion, or Z 2 x Z 2 
spin rotation symmetry [14, 24]. Bosonic SPT phases can be 
partially classified by group cohomology theory [25, 26]. Es¬ 
pecially, in ID (which we will study in this paper), SPT phases 
with onsite-symmetry are classified by 7^^(G, f/(l)), or the 
projective representations of the symmetry group G [27-29]. 

The ground-state degeneracy is a necessary condition for 
the existence of a ruled surface. However, in order to observe 
such a ruled surface, we need to add a local term that can lift 
the degeneracy. In other words, ground-state degeneracy that 
can be lifted by some local term will lead to ruled surfaces. 
It is then required such a local term exists. In case of the 
symmetry-breaking order, it corresponds to the local order pa¬ 
rameter. We show that in case of the SPT order, the emergence 
of a ruled surface only exists for system under open boundary 
condition (OBC), with the corresponding local term acting on 
the boundary that breaks the symmetry of the system, hence 
lifting the ground-state degeneracy. 

To show this, we study a ID model exhibiting Z 2 x Z 2 SPT 
order. We discuss in detail the effect of geometric locality 
and its relationship to the emergence of ruled surfaces. Since 
the degeneracy of the ground states is only exact in thermody¬ 
namic limit, in principle the ruled surface also requires such 
a limit. However, numerical results suggest that in practice, 
the ruled surface can already be observed for large but finite 
systems. This allows us to study the features of ruled surface 
based on finite-system calculations. 

One important difference between our results and Ref. [9] 
is that the boundary terms that lift the degeneracy are not asso¬ 
ciated with thermodynamic variables. It is essentially the ef¬ 
fect of geometric locality (i.e. boundary conditions do change 
the geometric locality of the system) that leads to a differ¬ 
ent geometry of the set of reduced density matrices. On the 
contrary, for a topological ordered system, no local terms can 
lift the topological degeneracy; therefore one cannot observe 
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Figure 1: ID system with A, B sublattice. The number of 
sites N is even. Green and blue dots represent sublattice A 
and B, respectively. Each small circle represents a unit cell, 
(a) is the ID system under open boundary condition and (b) 
is the system under periodic boundary condition. 

ruled surface on the geometry of local reduced density matri¬ 
ces. Our results hence lead to a deeper understanding of the 
physical meaning of the emergence of ruled surface. 

II. Z2 X Z2 SPT ORDER: THE ID CLUSTER STATE IN 
MAGNETIC FIELDS 

Consider a ID system with A, B sublattice under OBC, see 
Fig. 1(a). Assuming the number of sites N is even, the Hamil¬ 
tonian reads: 

N-l i=N 

i=2 i=l 

where X, Z are Pauli operators. 

The first term corresponds to the stabi¬ 

lizer generators of the IB cluster state (without boundary 
terms) [30, 31]. The second term corresponds to a longitudi¬ 
nal magnetic field. Notice that Hdu has a Z 2 x Z 2 symmetry 
generated by the following two operators (see e.g. [32-34]), 

Oi = ZiZsZ5...Z]sf-i ( 2 ) 

O 2 = Z 2 Z 4 ZQ...ZN (3) 

Obviously, [Oi, O 2 ] = 0 and [Oi.Hdu] = 0, [O 2 , Hdu] = 0. 
In fact, there is also a hidden continuous (7(1) symmetry in 
Hdu generated by 

When Bz < I, the ground state of above model has Z 2 x Z 2 
SPT order. To show this, we can transform Hdu into a fa¬ 
miliar form with only two-body interactions [35]. Consider 
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unitary operations Uab acting on each nearest A-B sites, as 
denoted in Fig.l. Each Uab transform Pauli operators as fol¬ 
lows: 


Xa Xa, Za ZaXb, (4) 

Xb Xb^ Zb XaZb- (5) 


In fact, Uab is nothing but the controlled-Z operation (i.e. 
CZ = diag(l, 1,1, —1)) in the Pauli X basis (i.e. ^(|0) ± 

|1))). That is, Uab = {Ha Z) Hb)CZab{Ha Hb), where 

H is the Hadamard transformation given hy HXH = Z. 
Under this transformation, we have 


and 


^ Xl^Zf (6) 

Xf_^Z^_Af (7) 

^ Zfx^^^ (8) 

^ Xfzf^^. (9) 


Thus we can recast the original Hamiltonian as following: 


N-2 

Hobc = Ji{^oi) {XiZi^i-\-ZiXi^i) 

i=2^even 
N-1 

+ J2(1— Q^) {XiZiZiXij^i). 

i=l^odd 


( 10 ) 


unchanged if the t/(l) symmetry is destroyed by anisotropic 
interactions. 

In later discussion, we will go back to the original clus¬ 
ter model. For the purpose of studying the reduced density 
matrix, we further add a transverse magnetic field. We will 
consider cases with both OBC and periodic boundary condi¬ 
tion (PBC). Thus under OBC (see Fig. 1(a)), the Hamiltonian 
reads: 


N-l 

Hobc = -^i(l + <^) ^ Xi-iZiXij^i 

i=2 

^ ^ (12) 

+ J 2 {l-a)Y,Zi-B^Y.^i 

= Ji(l + a)H° + J2(l - a)H° - B^H^, 
while under PBC (see Fig. 1(b)), the Hamiltonian is: 

AT+l 

Hpbc = Ml +a)Yl Xi-iZiXi+i 

i=2 

M ^ (13) 

+ J 2 (l — O') Zi — Bx Xi 

i=l i=l 

= Ji(l + a)Hf + J 2 (l - a)Hf - 

where the (A^ - 1 - 1, - 1 - 2) sites are identified with the (1, 2) 

sites, respectively. 


For the convenience of later calculation, we let Ji = ±1, 

J 2 = ±1, -1 < a < 1. If we further apply an unitary m. EFFECT OF LOCALITY AND THE EMERGENCE OF 
transformation on the A sublattice such that RULED SURFACE FOR SPT PHASE 


xf^zf, zf^xf, 

then the model (10) becomes the familiar XY model [36-39], 

N-2 

HoBC = Ml + oi) Y. {NXi+i + ZiZi+i) 

1=2,even 

N-l ^ ^ 

-\-J2{1 — Cl) {XiXij^i -\r ZiZij^i). 

i=l^odd 

The Z 2 X Z 2 symmetry becomes very clear in above Hamil¬ 
tonian: it is generated by the uniform X and Z operations. 
Since in the above model each unit cell contains two spins, 
the strong bonds may locate inter unit cells or intra unit 
cells, and consequently there are two different phases. When 
0 < a < 1, the ground state carries nontrivial Z 2 x Z 2 SPT or¬ 
der and has two-fold degenerate edge states (which carry pro¬ 
jective representation {/, X, Y, Z} of Z 2 x Z 2 group) on each 
boundary; on the contrary, when — 1 < a < 0, the model falls 
in a trivial symmetric phase without edge states; a = 0 cor¬ 
responds to the phase transition point. The previously men¬ 
tioned U{1) continuous symmetry is generated by Yi in 
Eq. (11). The U{1) symmetry is an accidental symmetry for 
the SPT order, namely, the properties of the two phase remains 


One necessary condition for the emergence of a ruled 
surface is the ground-state degeneracy. When B^ = 
0, (Ji, J 2 ) = (±l,zbl),0 < a < 1, the ground-state of 
HoBC is four-fold degenerate (if N -m 00 ) and the ground 
state of HPBC is unique. One may expect that a ruled surface 
will appear on the surface of the convex set Q{{Hj}) consist¬ 
ing of all the points given by 

MH?p), ^ p), ^ MMp)) (14) 

for any quantum state p, similar as the symmetry-breaking 
case as discussed in [9]. 

Unfortunately, it is not the case for the SPT order. This is 
because that the expectation value of the symmetry-breaking 
term iifs/X is mainly contributed from the bulk and the ruled 
surfaces which result from the edge states will become invis¬ 
ible. This is essentially the meaning of ‘topological’ for the 
SPT orders, in contrast to symmetry breaking orders. If one 
instead takes the expectation value of iifs, which is indeed dif¬ 
ferent for the degenerate ground states, the set 

(iV^ p), 1 tr(FOp),tr(ff3p)) , (15) 




4 


which is indeed convex, will be unbounded. Furthermore, if 
one only takes the expectation value on the boundary, i.e. to 
consider 

ltr(ffOp),tr((Xi +Xjv)p)) , (16) 

then this set is no longer convex (see Appendix A for more 
details). 

To overcome all these difficulties, we instead add the 
symmetry-breaking term on the boundary of the system. This 
is because that the degeneracy essentially comes from the 
edge spins. We therefore propose to use the following Hamil¬ 
tonian: 



Hobc = Ji(l + a)H? + J2(l - a)H° 
-B,{X^+Xn) 


(17) 


And for comparison purpose, we also modify the PBC 
Hamiltonian to be: 


HpBc = Ji{l + a)HP + J2{l-a)HP 

-B,{Xi+Xn) ’ 

For OBC, the convex set can be generated by the 

following expectation value with respect to the ground state: 


(XZX) = {H^)/{N - 2), 

(19) 

{Z) = {H^)/N, 

(20) 

(X) = {Xi + Xm)/2. 

(21) 

For PBC, the corresponding quantities are: 


{XZX) = {H[)/N, 

(22) 

{Z) = {HP)/N, 

(23) 

(X) = {X,+Xn)/2. 

(24) 

We show that there will then be emergence of ruled surfaces 
on the boundary of for OBC, and no ruled surfaces 

for PBC, which is illustrated in Fig. 2. 



IV. ALGORITHM AND PRECISION 

We study above models using two different approaches. We 
first study a small size system using exact diagonalization 
(ED) method, where there is no visible error. The numeri¬ 
cal result for = 12 is presented in Fig. 3, which already 
shows the signal of ruled surface under OBC. To go closer to 
the thermodynamic limit, we use matrix product state (MPS) 
as a variational ansatz and approach the ground state using 
Time-Evolving Block Decimation (TEBD) method [40-42], 
whose accuracy is mainly limited by Trotter error and finite 
dimension of the underlying MPS. In practical calculations, 
the ruled surface can be distinguished from a non-ruled one 
by the oscillating scenario in the convex set which arises due 
to the ground-state degeneracy, as further discussed in the next 



Figure 2: Schematic figure for the convex set under OBC and 
PBC. Red dots in (a) mark the phase transition point. The 
green surface in (a) is the ruled surface, which is absent in 
(b). For clarity, only the upper half (i.e. {X) >0) is shown. 
For more details, see the analysis in Sec. IV. 


paragraph. As shown in Fig. 4, the green oscillating line indi¬ 
cating the ruled surface indeed only exists under OBC. 

In the thermodynamic limit under OBC, when the system 
is in the SPT phase, the ground state space is spanned by 4 
degenerate states. For a large finite system, these four states 
are nearly degenerate. Thus the state given by TEBD method 
would be a superposition of these four states because of the 
limit of the numerical accuracy. This explains the vibrational 
property of the ruled surface in Fig. 4 (a). As far as we are 
mainly concerned with the extent of the ruled surface, it is 
safe to replace the original vibrating curve with its upper hull. 
Similar vibration was also observed when external magnetic 
field Bx is small enough, e.g. 10“^, which is not shown in the 
figure. 

A thorough investigation of the numeric errors would be 
both lengthy and unnecessary. Here we perform a qualitative 
analysis to show how such errors interestingly lead to the pos¬ 
sibility to obtain the ruled surface in a large but finite system, 
while such a surface should only exist in the thermodynamic 
limit. 

Due the the limit of numerical accuracy, the curve com¬ 
puted for Bx = 0 should be more properly understood as 
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Figure 3: Convex set for ID cluster model with system size 
= 12 . (a) represents convex set under OBC while (b) is 
for PBC. For both cases, the field is added only on two 
boundary sites (one for each boundary). The result is 
obtained using ED method. Bx =0 is represented by the 
green line in both cases. A dramatic difference between the 
two cases is that a small region reminiscent of ruled surface 
exists in (a) but is absent in (b). 


the curve corresponding io Bx = e, where e is a very small 
number. Denote by fN,B^ the curve for a system of size N 
and the external magnetic field Bx, or the upper hull of it in 
case of vibration. Thus our goal is to estimate the difference 
between /oo,o (theoretical boundary for the ruled surface in 
the thermodynamic limit) and f]sf,e (curve observed in a fi¬ 
nite system with numerical errors). Same as finite system, we 
can take /oo,o ~ foo,e- Since foo,e = limAr^oo fN,e, for any 
(5 > 0, there exists N{6), such that d(/oo,e,/Ar,e) < S for 
any N > N{6), where d can be taken, for example, to be the 
Hausdorff distance between curves. Thus the difference be¬ 
tween /oo,o ~ /oo,e and fN,o ~ fN,€ can be arbitrarily small 
for N large enough. In practice, the convergence is fast such 
that when N = 60, the observed ruled surface precisely rep¬ 
resents the ruled surface in the thermodynamic limit. 

For a large finite system under PBC, there seems to be a 
small ruled surface at the phase transition point, shown in 
Fig. 4 (b). With increasing bond dimension D, this small area 
shrinks and eventually vanish in the infinite D limit, shown in 


Figure 4: Convex set for ID cluster model with system size 
= 60. (a) is for OBC while (b) is for PBC. The green 
vibrational curve is the ruled surface. We use TEBD method 
with internal dimension D = 40. Due to numerical error, we 
can witness two ruled surface with a large but finite system 
under OBC. 


Fig. 5. Thus under PBC, there is no ruled surface. 

Notice that under both OBC and PBC the upper plane is 
flat. This is because the normal direction of the corresponding 
supporting hyperplane is (0, 0,1), which corresponds to Ji = 
J 2 = 0 for Hqbc HpBC' Hqbc Hpbc then both 
become —Bx{Xi Xn), which only acts nontrivially on the 
boundary, hence are largely degenerate. On the contrary, each 
line inside the ruled surface (only under OBC) corresponds to 
a (finite) four-fold degeneracy, which is a non-trivial signal of 
the SPT order. 


V. CONCLUSION AND DISCUSSION 

We study geometry of reduced density matrices for SPT or¬ 
der. Our focus is on the emergence of ruled surface on the 
boundary of the convex set The ground-state de¬ 

generacy is a necessary condition for the existence of those 
ruled surfaces, yet not sufficient. 

Compared to the ruled surfaces associated with symmetry¬ 
breaking order as discussed in [9], there is an essential dif- 
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(a) 



Figure 5: Finite D scaling for phase transition point under 
PBC. In both (a) and (b), different colors represent different 
bond dimension D where yellow, blue and green is for 
= 40, 60, 80, respectively. Filled square represents 
Bx = I X while filled circle represents Bx = 0. With 
increasing D, the small ruled surface in Fig. 4 (b) near phase 
transition point shrinks. Eventually, this small ruled surface 
will vanish in the infinite D limit. 


ference for the SPT order. Since there is no local order pa¬ 
rameter for SPT order, the ruled surface only exists for the 
open boundary condition with symmetry-breaking term acting 
on the boundary. This term is not a thermodynamic variable. 
Therefore the emergence of ruled surface for SPT order is an 
effect of geometric locality of the system. 

In principle, ruled surface only exists in the thermodynamic 
limit. However, we have shown that in practice, finite-size cal¬ 
culation suffices to reveal this phenomenon, due to inevitable 
computational precision uncertainty. This allows us to deal 
with the calculations using finite systems. 

We hope our discussion leads to further understanding of 
the geometry of reduced density matrices, the effect of geo¬ 
metric locality, and SPT order. 
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Appendix A: Vanishing of ruled surface and nonconvex set 

In this Appendix, we explain in detail why we choose the 
Hamiltonian (17) instead of (12). 

If we choose the Hamiltonian (12) and plot the convex set 
(14), the ruled surface will vanish in the thermodynamic limit. 
The reason is that the degeneracy of the ground states are ow¬ 
ing to the edge states, while the expectation value of the term 
^tr(i73p) mainly comes from the bulk. If ^ oo, the 
boundary effect (together with the ruled surface) will disap¬ 
pear due to the normalization factor ^ (See Eig. 6). 



Eigure 6: Convex set with a small ruled surface, which is 
indicated by green lines. The result is obtained by ED 
method with AT = 12. {XZX) = ), 

On the other hand, if we plot the set (16) to avoid the 
vanishing factor what we obtain is not a convex set (see 
Eig. 7). This is because we cannot use {Hi ^ H 2 ^Xi X^} 
to construct the Hamiltonian (12). 














7 



Figure 7: Non-convex set. Result is obtained by ED method 
withiV = 12. {XZX) = ), {Z) = 

{X) = {Xj+Xm). 
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